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The fundamental parameters approach to 
line profile fitting uses physically based 
models to generate the line profile shapes. 
Fundamental parameters profile fitting 
(FPPF) has been used to synthesize and fit 
data from both parallel beam and divergent 
beam diffractometers. The refined parame- 
ters are determined by the diflfractometer 
configuration. In a divergent beam diflfrac- 
tometer these include the angular aperture 
of the divergence sht, the width and axial 
length of the receiving slit, the angular 
apertures of the axial Soller slits, the 
length and projected width of the x-ray 
source, the absorption coefficient and axial 
length of the sample. In a parallel beam 
system the principal parameters are the 
angular aperture of the equatorial 
analyser/Soller slits and the angular aper- 
tures of the axial Soller slits. The presence 
of a monochromator in the beam path is 
normally accommodated by modifying the 



wavelength spectrum and/or by changing 
one or more of the axial divergence 
parameters. Flat analyzer crystals have 
been incorporated into FPPF as a 
Lorentzian shaped angular acceptance 
function. One of the intrinsic benefits of 
the fundamental parameters approach is its 
adaptability any laboratory diflfractometer. 
Good fits can normally be obtained over 
the whole 20 range without refinement 
using the known properties oflthe diflfrac- 
tometer, such as the sUt sizes and diflfrac- 
tometer radius, and emission profile. 
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1. Introduction 

The fundamental parameters approach to line profile 
fitting uses physically based models to generate the line 
profile shapes. The instrument profile shape K(2Q) is 
first synthesised by convoluting together the geometri- 
cal instrument function J(2Q) with the wavelength pro- 
file WilQ) at the Bragg angle 26^ of the peak, 

K{2Q) = jWi20-2<p)Ji2<p)d2(p = Wi20) ® Ji20) (1) 

where the function J{20) itself is a convolution of the 
various instrument aberration functions associated with 
the diffractometer. 



k„Ji2O)=J,i2O)0J2i2e)0,,,0Ji2e) ®Jn(2^). 

(2) 

Diffraction broadening is incorporated into the profile 
fiinction 1(20) by convoluting the broadening fiinction 
B(2$) into the instrument profile fiinction as shown Fig. 
1, 



i(2e)=K(2e)®B(2e). 



(3) 



This technique of profile synthesis was first introduced 
50 years ago by Alexander [1], but has only been imple- 
mented as a standard fitting procedure during the last 
ten years [2,3,4]. More recently, freeware and commer- 



Deceased. 
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Fig. 1. Convolution of the geometric instrument profile J(20)^ the wavelength profile W(20) and the specimen 
broadening function B{29) to produce the measued convoluted function I{20), 



cial software packages [5,6,7] have become available 
for fundamental parameters profile fitting (FPPF) 
either for use as single line profile fitting, lattice param- 
eter refinement or for Rietveld analysis. 

FPPF has been used to synthesise and fit data from 
both parallel beam and divergent beam diffractometers. 
The refined parameters are determined by the diffrac- 
tometer configuration. In a divergent beam diffrac- 
tometer these include the angular aperture of the diver- 
gence slit, the width and axial length of the receiving 
slit, the angular apertures of the axial Soller slits, the 
length and projected width of the x-ray source, the 
absorption coefficient and axial length of the sample. In 
a parallel beam system the principal parameters are the 
angular aperture of the equatorial analyser/Soller slits 
and the angular apertures of the axial Soller slits. The 
presence of a monochromator in the beam path is nor- 
mally accommodated by modifying the wavelength 
spectrum and/or by changing one or more of the axial 
divergence parameters. Flat analyser crystals have been 
incorporated into FPPF as a Lorentzian shaped angular 
acceptance function. 

One of the intrinsic benefits of the fundamental 
parameters approach is its adaptability to any laborato- 
ry diffractometer. Good fits can normally be obtained 
over the whole 20 range without refinement using the 
known properties of the diffractometer, such as the slit 
sizes and diffractometer radius, and the emission pro- 
file. Fine tuning is sometimes necessary to accommo- 
date a monochromator or to compensate for the fact 
that certain aberrations are not completely independent 
[8]. Under these conditions some of the instrument 
parameters need to be refined, but the refined values 
normally are within ±10% of the actual values. 
Correlation between refined instrument parameters can 



occur when fitting to data over a restricted 20 range. 
Such correlation occurs between the axial divergence 
parameters and absorption as both of these aberrations 
can produce similar forms of asymmetric profiles 
between 20= 50° and 100° in diverging beam diffrac- 
tometers. Correlation is minimised by using data with a 
large 20 range so that the unique angular dependence of 
individual aberrations becomes evident. When a set of 
instrument profiles cannot be fitted by FPPF, this is 
usually an indication that either the model used is 
invalid (eg. incorrectly chosen slit value), the instru- 
ment is mis-aligned, there is an overlapping impurity 
line or, the specimen is generating crystallite size 
broadening or is inhomogeneously strained. 

FPPF was designed originally as a tool for analysing 
diffraction line broadening. Fitting is done by convolu- 
tion and corrections for instrument broadening and 
peak shift are intrinsic to the refinement. When an 
instrument is well characterised, line broadening can be 
analysed without a reference specimen. Moreover, 
when a reference standard is used, which has different 
properties from the specimen with line broadening, 
some compensation can be made for these differences. 
For example, when LaB^ SRM 660a (//powder ^ 500 cm"^) 
is used as a reference, compensation can be made for 
differences in the absorption of the sample and LaB^. In 
the latest version of the commercial software package 
(TOPAS)^, the concept of fundamental parameters has 
been extended so that any user defined profile that 



Certain commercial equipment, instruments, or materials are iden- 
tified in this paper to foster understanding. Such identification does 
not imply recommendation or endorsement by the National Institute 
of Standards and Technology, nor does it imply that the materials or 
equipment identified are necessarily the best available for the pur- 
pose. 
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accurately describes the physical broadening can be 
readily convoluted into the refinement. 

In this paper we will discuss the physical origin of 
the instrumental profile shapes for various laboratory 
diffractometer configurations including both divergent 
beam and parallel beam instruments. This will include 
a description of the geometrical aberrations as well as 
discussion on the nature of the wavelength distribution 
and the influence of monochromators on this distribu- 
tion. Some discussion is also presented to demonstrate 
how the FPPF may be fitted to experimental data from 
a material with a low attenuation coefficient. 



2. Basic Objectives of the FPPF 
Technique 

One of the basic objectives of the FPPF technique is 
to be able to fit any powder diffraction profile using a 
physically based model to describe both the instrument 
profile and any diffraction broadening generated by the 
specimen. In principle, therefore, the technique should 
be adaptable to any powder diffractometer and fit pro- 
files of widely differing shapes, such as those in Fig. 2, 
by simply modifying the physical parameters of the dif- 
fractometer used to describe the profile. 

Although most applications of FPPF have focussed 
on the conventional diffractometer it has also been 
utilised for analysing neutron diffraction data [9] and 



synchrotron data [10,11]. In the TOPAS implementa- 
tion of FPPF there are a wide variety of possible aber- 
ration functions available within the package and these 
can put together to suit a particular diffractometer 
design and in terms of parameters that are relevant to 
the instrument. One of the most important achieve- 
ments of the FPPF technique for practical users is speed 
of calculation. Accurate multiple convolution calcula- 
tions over large 20 ranges can be very time intensive 
and it is of central importance to minimise this time to 
enable Rietveld refinement to be completed in "sec- 
onds" without loss of accuracy within the profile func- 
tion synthesis procedures. Various procedures have 
been implemented, some of which are described by 
Cheary and Coelho [2,3,4], but one of the most impor- 
tant has been to code the time intensive calculations at 
an assembler code level taking steps to optimise the use 
of the various registers within the PC chip. 



3. Laboratory Diffractometer 
Configurations and Their Geometrical 
Aberrations 

Up until the mid-1990s most of the laboratory pow- 
der diffractometers in use were divergent beam instru- 
ments with a narrow receiving slit, diffracted beam 
monochromator and a simple proportional/scintillation 
counter detector as shown in Fig. 3a. Over the past 10 
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Fig. 2. Comparison of profiles produced by the reference specimen LaB^ (SRM 660a) at sim- 
ilar 29 angles on two different diffractometers, (a) on the high resolution powder diffractome- 
ter BM 16 at the ESRF, Grenoble (A =0.35 A) and, (b) on a conventional ^20 divergent beam 
diffractometer (A = 1.541 A) with a radius of 217 mm and commonly used sht sizes (ie. 1° 
divergence sHt, 0.2 mm receiving slit and 2° SoUer slits in the incident beam). 
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Fig. 3. Two configurations of diverging beam diffractometer showing the princi- 
pal optical components, (a) with a diffracted beam monochromator and, (b) with 
an incident beam monochromator. 



years the number of diffractometer options available 
from manufacturers has increased and users operate 
with a wider range of x-ray optical designs. The types 
of geometrical aberrations encountered is somewhat 
broader than those discussed in the classic work by 
Wilson [12] as will be discussed below. 

3.1 Divergent Beam Diffractometers — Symmetric 
Diffraction 

The most widely used laboratory diffractometer in 
use today is still the divergent beam diffractometer with 
either a bent graphite monochromator in the diffracted 
beam (see Fig. 3a) or, a ground and bent asymmetrical- 
ly cut germanium monochromator in the incident beam 
(see Fig. 3b). Both of these configurations possess a 
similar array of geometrical aberrations. The major dif- 
ference between them is the wavelength distribution 
which normally consists of both the Kai and Ka2 com- 
ponents of the K spectrum in the graphite monochro- 
mator case and only the Ka^ with the Ge monochroma- 
tor [13]. Further discussion of the wavelength distribu- 
tion and the effects of monochromators is given later. 

The principal geometric aberrations contributing to 
profiles from the above diffractometers are, 
(i) the finite width of the x-ray source, 



(ii) a divergent incident beam on to a flat specimen 

(flat specimen error), 
(iii) the finite width of the receiving slit, 
(iv) the beam penetration into the specimen (specimen 

transparency), 
(v) the deviation of the beam from the equatorial plane 

(axial divergence). 
These aberrations all produce some degree of line 
broadening and, in the case of flat specimen error, spec- 
imen transparency and axial divergence, some asym- 
metry is also introduced. Zero 2^ and specimen surface 
displacement errors may also be present in a diffrac- 
tometer, but these only affect the 20 position of a pro- 
file and not its shape. In both configurations the mono- 
chromators not only determine the wavelength distribu- 
tion, but they also act to reduce the axial divergence 
and it is often considered unnecessary to include Soller 
slits between the sample and the monochromator. 

3.2 Divergent Beam Diffractometers — 
Asymmetric Diffraction 

Divergent beam diffractometers used under symmet- 
ric conditions only measure diffraction from planes 
parallel to the specimen surface. To measure diffraction 
from planes angled relative to the specimen surface it is 
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necessary to operate under asymmetric conditions as 
illustrated in Fig. 4a. The problem with operating in 
this mode on most commercial diffractometers is that 
the receiving slit is no longer at the focus of the dif- 
fracted beam and profiles are broadened by "defo- 
cussing". The amount of defocussing is determined by 
the angle of divergence of the incident beam and dis- 
tance of the focus from the receiving slit. 

Defocussing also occurs in diffractometer configura- 
tions where, 

• the sample is oscillated ±50) about the diffractometer 
axis so that the angle of incidence on to the specimen 
varies between -^ 5(0 and - 5(0. This oscillation 
moves the focus of the diffracted beam continuously 
back and forth in front of and behind the receiving 
slit, 

• the receiving slit is replaced by a position sensitive 
detector (PSD). The only position on the detector that 
is normally in focus is its centre (see Fig. 4b); all dif- 
fracted beams entering the detector at off-centre posi- 
tions are defocussed. In PSD systems, the aberrations 
contributing to a profile include all the standard dif- 
fractometer aberrations, except the aberrations for- 
merly due to the receiving slit are replaced by defo- 
cussing, the discharge resolution of the detector, and 
parallax error [14,15]. In a scanning PSD diffrac- 
tometer, the recorded profile shape is an average of 
all the profile shapes across the active window 
length. 



3.3 Parallel Beam Diffractometers 

There are two common forms of the parallel beam 
powder diffractometer which are illustrated in Fig. 5. 
These are based on using either analyser slits (other- 
wise referred to as equatorial Soller slits) or a flat Ge/Si 
analyser crystal as the angular discriminator of the dif- 
fracted beam. Amongst laboratory diffractometers the 
analyser slit set-up is the most widely used form as it 
offers more intensity but poorer resolution than the 
analyser crystal set-up. 

Parallel beam diffractometers have emerged as one 
of the most popular forms of laboratory diffractometer 
over the past ten years and now constitute more than 
30 % of the new diffractometer purchases. There are 
fewer geometric aberrations contributing to the profile 
and systematic errors arising from specimen displace- 
ment, specimen transparency and surface roughness are 
not significant. There are two geometric aberrations 
contributing to a parallel diffractometer, 
(i) the angular acceptance function of the analyser 

foils or analysing crystal, 
(ii) deviation of the beam from the equatorial plane (ie. 

axial divergence). 
In most laboratory diffractometers, the parallel beam is 
produced by using a parabolic graded multilayer mirror 
with the line x-ray source positioned at the focus of the 
mirror [16]. Although the beam may be parallel in the 
equatorial plane, it will not be parallel in axial plane 
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Fig. 4. (a) Diverging beam diffractometer operated in asymmetric mode with defocussed dif- 
fracted beam at the receiving slit, (b) diffractometer with PSD replacing the receiving slit. 
Diffracted beam is in focus at the centre of the detector and defocussed in off-centre positions. 
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Fig. 5. Two configurations of parallel beam diffractometer, (a) using analyser slits 
in the diffracted beam and, (b) using a flat analyser crystal. 



and axial divergence can be expected in both the inci- 
dent and diffracted beams. Low angle profiles will 
therefore be asymmetric although not to the same 
extent as diverging beam instruments. 



4. The Instrument Aberrations 

The geometric instrument aberrations tend to deter- 
mine the shape of a diffractometer profile at low 20 
angles (ie., 20 < 50^). At high 20 angles (20 > 100°), 
the profile conforms primarily to the shape of the wave- 
length distribution in the beam. With the exception of 
the aberrations associated with "receiving system" of 
the diffractometer and the x-ray source, all of the geo- 
metric instrument aberration profiles vary with 20. In 
the following sections the shapes of the major instru- 
ment aberrations used in FPPF analysis to describe the 
various laboratory diffractometer configurations are 
discussed for conditions that are typical of those 



encountered in practice. The aberration functions gen- 
erated by mis-setting a diverging beam diffractometer 
or using it under asymmetric conditions are also dis- 
cussed. Most of the results quoted here are for a diffrac- 
tometer radius R = 2\5 mm. 

The convention adopted here for describing the 
angular variables is that 20 refers to the continuously 
variable angle measured on the diffractometer whereas 
2^ or 20^ refers to the Bragg angle of the diffraction 
line. The angle £ refers to the difference between the 
measured angle 20 and 20, 



£ = 2(l)-20 
4.1 Finite X-Ray Source Width 



(4) 



The profile shape of this aberration is generally 
expressed as an impulse function of width A20^ as 
shown Fig. 6a. Although the choice of an impulse func- 
tion may be not be strictly valid for describing the x-ray 



(a) 



AZOx = Wx/R rad 
< — ► 



(b) 



H^ £ 



Fig. 6. Aberration profile models, (a) simple model for a source of projected width w^ and a diffractometer radius R, (b) 
model with "tube tails" containing additional parameters/^ /tai/Zniax ^^^ angular widths Z^ and Z2 from the central maximum. 
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source aberration, the exact shape used is not critical 
when a long fine focus tube (target width -0.4 mm) is 
installed on the diffractometer. At a take-off angle of 6° 
this appears as a projected width w^ ~ 0.04 mm and the 
aberration profile has a width A20^ -0.01° and does not 
contribute significantly to the overall width of the 
instrument profile. 

In broad focus tubes, the target width -2 mm and 
w^ - 0.2 mm at 6° take-off so that the aberration profile 
width A20^ - 0.056°. At this level the source width 
makes a much bigger contribution to the overall width 
of the instrument profile and a more accurate form for 
the aberration profile shape is necessary. A good 
approximation under these circumstances is a Gaussian 
shape rather than an impulse function. In diffractome- 
ters with curved crystal incident beam monochromators 
the source width can also have a greater contribution 
because of the magnification effect introduced by the 
monochromator. This occurs with asymmetrically-cut 
Johansson incident beam monochromators where the 
source-crystal and crystal-focal point distances are typ- 
ically -120 mm and -230 mm, respectively. A fine 
focus tube with a projected width of 0.04 mm is then 
effectively magnified to -0.08 mm. Under these condi- 
tions the effective source width can be trimmed down 
by reducing the width of the focal line slit. 

For accurate line profile analysis it is also necessary 
to modify the simple impulse model even with long 
fine focus tubes. Bergmann [ 1 7] has shown that most of 
the anode surface in an x-ray tube produces x rays 
albeit at a much lower intensity than the focal line on 
the anode. This is illustrated in Fig. 7 which shows the 
intensity recorded by scanning with a 50 |Lim slit across 
the image of a x-ray source formed through a 10 |Lim 
pinhole in platinum. A better approximation to the aber- 
ration function is a sharp impulse function superim- 
posed on a broad impulse function to represent the so 
called "tube tails". This is illustrated in Fig. 6b. The 
parameters introduced to describe the "tube tails" are 
the extents of the high and low angle tails, Zj and Z2, 
and the intensity of the tail f relative to the intensity at 
the tube focus. In most instances the intensity of the 
tails is -0.1 % of the peak intensity and is only signifi- 
cant when analysing intense lines. The tails themselves 
are not necessarily symmetric with respect to the tube 
focus and can extend over a 20 range up to 0.6°. 

4.2 X-Ray Receiving System Models 

In diverging beam diffractometers the receiving slit 
is placed at the focus of the diffracted beam and for per- 
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Fig. 7. Intensity scan with 50 |ani wide slit of an image formed 
through a 10 frm pinhole in platinum of the 0.4 mm wide long fine 
focus in a Cu anode x-ray tube set at 40 kV, 40 mA. 



feet focussing should have an infinitely small width. 
Owing to the many aberrations present, focussing is 
never perfect and the count rate incident on the receiv- 
ing slit tends to increase with increasing slit width, but 
at the expense of resolution. In parallel beam diffrac- 
tometers the receiving system is based on using either 
the Hart-Parrish system of analyser slits [18] or, a flat 
analyser system as illustrated in Fig. 5 earlier. In many 
glancing incidence diffractometers the receiving sys- 
tem consists of analyser slits and an analyser crystal in 
the diffracted beam. Although the aberration functions 
associated with the various receiving systems for paral- 
lel beam and divergent beam diffractometers all pos- 
sess different shapes, they all possess the common 
property of being independent of 2^. 

4.2.1 Receiving Slit in a Diverging Beam 
Diffractometer 

Most commercial diffractometers have a selection of 
receiving slits ranging in width from 0.05 mm up to 0.3 
mm although occasionally larger slit sizes up to 0.6 mm 
are used to measure integrated intensity rapidly. The 
aberration function for a perfectly aligned receiving slit 
is an impulse function of width A20, given by 



A2e =^rad 
' R 



(5) 



where w^ is the width of the receiving slit. The angular 
width A2^, subtended by the receiving slit relative to 
the diffractometer axis is therefore normally between 
0.013° (0.05 mm) and 0.08° (0.3 mm). When the slit 
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size is larger than 0. 15 mm, the receiving slit aberration 
is often the dominant aberration in a diffractometer 
over the angular range 20= 15"^ to 60"^. 

4.2.2 Parrish-Hart Analyser Slits 

Analyser slits act as an angular filter in the diffract- 
ed beam. The aberration function or transmission fiinc- 
tion for these slits is a triangle fiinction, as shown in 
Fig. 8, in which the base width A20, is given by the 
angular aperture A of the slits. 

In the original Parrish-Hart diffractometer on Station 
2.3 at the Daresbury synchrotron, the analyser slits 
were 360 mm with a spacing of 0.2 mm between adja- 
cent foils giving an angular 20 aperture A~ 0.06°. In 
laboratory diffractometers the angular aperture A is typ- 
ically -0.1°. A problem often encountered with 
analyser slits is specular x-ray reflection from the 
analyser foils [19]. Weak satellite peaks appear on both 
the high angle and low angle profile tails but not neces- 
sarily of the same intensity as shown in Fig. 9a. This 



effect has been incorporated into the aberration profile 
by adding two Voigt functions of unequal intensity ,one 
on each side of the triangular aberration function, to 
represent the satellite reflections [11]. The parameters 
of the satellite peaks can then be determined by fitting 
profiles from a reference material such as the NIST ref- 
erence standard LaB^, SRM 660a. An alternative 
approach to determining the aberration profile of an 
analyser slit system, without the effects of the wave- 
length profile distorting the result, is to simply carry 
out a 2^ scan across the incident beam. Provided the 
axial divergence of the incident beam is kept small, by 
including axial Soller slits, and the equatorial diver- 
gence is negligible then the incident beam scan will 
have exactly the same shape as the aberration profile of 
the analyser slits. 

4.2.3 Analyser Crystals 

The inclusion of an analyser crystal in the diffracted 
beam of a parallel beam diffractometer gives high res- 
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Fig. 8. Triangle shaped aberration function for a set of analyser slits with an angu- 
lar aperture A where A/2 = spacing between the foil/length of the foils. 
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Fig. 9. (a) Reflection satellite peaks from analyser slit recorded using the 310 line from NIST standard material LaB6 SRM 660a using the 
diffractometer on Station 2.3 at Daresbury synchrotron, (b) 20 scan across a 0.1 mmx 0.1 mm incident beam using a Gel 11 analyser crystal 
on beamline BM16 at the ESRF, Grenoble. 
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olution diffraction patterns with a low background, but 
the intensity is invariably less than the Parrish-Hart 
configuration. The aberration profile introduced by the 
analyser crystal is generally very narrow and can be 
determined by measuring the rocking curve of the crys- 
tal. For a perfect analyser crystal the aberration profile 
will be determined by the Darwin profile of the 
analyser crystal. In practice, however, the aberration 
profile will be broadened by the mosaic structure of the 
crystal, any stresses in the crystal and any waviness or 
curvature of the crystal surface [20]. As a consequence 
the aberration profile can be dependent on the size of 
the beam incident on the crystal. A first approximation 
to the shape of the aberration profile of an in-situ 
analyser can be obtained from a 20 scan of the 
analyser/detector using a very fine incident beam as 
shown in Fig. 9b. Although the profile recorded in this 
way also has the wavelength distribution folded into it, 
the result does at least give an indication of the shape 
and upper limit of the FWHM of the aberration profile. 

4.3 Flat Specimen Error 

The basic optics of the focussing powder diffrac- 
tometer set up for symmetric diffraction is illustrated in 
Fig. 10. The x rays are incident at an angle ^ on an ideal 
polycrystalline specimen with a surface radius of cur- 
vature p. For diffraction from a particular hkl plane the 
common property of all the diffracted rays from the 
specimen is that they all deviate through the same angle 
20. By simple geometry it can be shown that all the dif- 
fracted rays converge to a focus on a circle which has 
the same curvature as the specimen surface. The focus 
of the diffracted rays defines the position of the receiv- 

Curved Polycrystalline 
Specimen 



X-ray 
Source 




Receiving 
Slit 



ing slit. In commercial diffractometers the specimen is 
invariably flat and the diffracted beam no longer 
focusses perfectly. Good focussing characteristics, 
however, can be maintained with appropriately chosen 
slits to limit the equatorial divergence and a diffrac- 
tometer radius R sufficiently large to reduce defo- 
cussing errors to an acceptable level without losing too 
much diffracted intensity. 

For an incident beam, with an equatorial divergence 
a centred on the diffractometer axis, the aberration pro- 
file J^^(20) is asymmetric and exists only for the region 
£= 0. X rays diffracted from the centre of the specimen 
are detected at 20 = 2^ where as x rays diffracted off- 
centre are detected at 20 < 2^ as shown in Fig. 1 1 . 




Fig. 11. Diffraction from a flat plate showing the relationship 
between the measured angle 2(j> on the diffractometer and the diffrac- 
tion angle 20 for a ray at the outer limit of a beam of divergence a. 



The relationship between the difference £=2^-20 
and the distance q from the diffractometer axis at which 
the ray is diffracted is. 



e = -\^\ sin2erad 



(6) 



assuming small angles of divergence for the incident 
beam (eg. a < 2°). Thus, when a specimen is illuminat- 
ed over a length Q then the x rays diffracted from either 
extremity of the beam (ie., at ^ = ^QI2) defines the lim- 
iting value of £m of the aberradon function JpsC^) which 
exists within the range > £ > £m where 



2R 



sin 20 rad 



(7) 



This relation can also be expressed in terms of the 
equatorial divergence of the beam a starting from the 
relation. 



,=^f. 



1 



1 



2 lsin(e-a/2 sin(e+a/2 



(8) 



Fig. 10. Focussing in a synnnetric powder diffractometer. 
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When 20> 10°, Eq. (8) can be approximated as 
aR 



Q- 



sin^ 



and Eyi becomes 



£., = — cot 6^ rad. 

M 2 



(8a) 



(9) 



When the incident beam is centred on the diffractome- 
ter axis the normahsed equation for the aberration func- 
tion Jfs(^) f'or fl^t specimen error is, 



Jps(£): 



1 



2Jee, 



for > £ > ^t^. 



(10) 



Modem commercial diffractometers operate with either 
a fixed divergence a or a fixed illumination length Q. 
In the fixed a mode the aberration function is broad at 
low 20 and £m has a cot dependence whereas in fixed 
Q mode the breadth rises from zero at 2^= up to a 
maximum at 2^ = 90°. The extent of the changes in the 
aberration function Jsf(^) for each mode of operation 
using typical operating values for a and Q are shown in 
Fig- 12. 

With a fixed angle of divergence a = 1 °, the effects 
of flat specimen error in commercial diffractometers 
are discemable as an increase in both the asymmetry 
and breadth below 20^40° [21]. Under conditions of 
constant a the beam size increases with decreasing 20 
and eventually the beam will cover the whole speci- 
men. The angle 20^^^ at which this occurs is given by. 



smiO^) = aR/L 



(11) 



where L^^ is the length of the specimen. Under these cir- 
cumstances the value of £m is given by Eq. (7) with 



Q = Z,p. Below 20^1^ the aberration function remains the 
same as the beam extends beyond the specimen. 

In diffractometers with a fixed illuminated length, 
the effects of flat specimen error are generally smaller 
at low 20 values and increase with increasing 20. With 
an illuminated length of 20 mm on the specimen, flat 
specimen error is clearly discemable at 2^> 20°. It is 
not always possible to collect data from fixed Q mode 
diffractometer over a large 2^ range (eg., up to 150° 20 
as the a angle required at large 2^ is larger than the dif- 
fractometer can accommodate. For example, to main- 
tain a fixed beam length of 20 mm on a specimen over 
the range 2^ = 0° to 90°, the angle of divergence a will 
need to increase from 0° up to -4°. A divergence angle 
a = 4° is close to the maximum value at which most 
diffractometers can operate when a pyrolytic graphite 
monochromator is installed in the diffracted beam. At 
angles of a greater than 4° the diffracted beam may 
extend beyond the graphite crystal and not be diffract- 
ed into the detector. In practice, this is not usually a 
problem as fixed Q mode diffractometry is normally 
used for analysing materials such as clays with very 
low angle diffraction lines where the lines of interest 
start at 20 -3°. 

4.4 Specimen Transparency 

Specimen transparency produces asymmetry and 
broadening of the instrumental profile function. For 
perfect focussing all the diffraction should occur on the 
focussing circle, but when the beam penetrates the sur- 
face, diffraction will occur over a range of depths with- 
in the specimen. The aberration function J^ (e) for an 
infinitely thick specimen is given by 



-/„(£)= 



exp{e/d) 



e<0 



(12) 




19.85 



19.90 



19.95 



20.00 



20.05 



a= I 



Q =20 mm 




■-■ St^ 



59.85 



59.90 



59.95 



60.00 



60.05 



Fig. 12. Comparison of flat specimen aberration functions at 20= 20° and 20=60° for a diffractometer with either a fixed diver- 
gence angle a= 1° or fixed illuminafion length Q = 20 mm (7? = 215 mm). All the aberration functions shown were convoluted with 
a very narrow Lorentzian emission profile with a FWHM = 0.001 mA to overcome infinities at e= 0. Each plot is normalised to the 
same I^^^ and both plots cover a range of 0.2° 29. 
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where 5 = {lljxR) sin2^rad and jx is the linear attenua- 
tion coefficient. In low absorbing specimens which 
cannot be considered to be infinitely thick, the angular 
variable £ has a lower limit e^^ so that in Eq. (12), 



0>f;>£ 



and 



7T_ 
R 



cos ^rad 



(12a) 



where T is the specimen thickness. As a consequence 
the unit area normalising constant in (12) also changes 
and J^(£) becomes, 

exp(£ / 5) 



M^) = 



5[l-exp(£^/5)]' 



(12b) 



The asymmetry and broadening from specimen 
transparency is greatest for low absorption materials 
and is clearly evident when the linear attenuation coef- 
ficient jU < 50 cm~\ The contribution of specimen trans- 
parency is greatest at 2^-90° and at this angle the 
aberration profile has a FWHM - 0.03° 2^ when ^u - 50 
cm"\ but this drops to --0.005° 20 when jx - 200 cm"\ 
Fig. 13 shows the shapes of aberration profiles for an 
infinitely thick specimen and how they are affected by 
both jU and 2ft Specimen transparency effects are quite 
strong in polymeric materials with very low attenuation 
coefficients (jU ^ 30 cm"^ or less). They can also show 
up in loosely bound powders of low atomic number 
materials (e.g., MgO or Si) where the porosity ^50 % 
or less so that a material with a jU = 100 cm"^ is reduced 
to a powder with a jU = 50 cm"\ 

4.5 Diffractometer Defocussing 

Defocussing results in broadened diffiaction lines 
and occurs when the receiving slit is not positioned at 
the focus of the diffiacted beam. The most common 
causes of defocussing are, 

• mis-setting the incident beam angle O) so that it is no 
longer at the symmetric condition 6> = ft as in asym- 
metric diffiaction or. 



• wrongly positioning the receiving slit so that the dis- 
tance of the slit from the sample is larger or smaller 
than the nominal radius R of the diffractometer. 
For both of these conditions the focus of the diffiacted 
beam will be either in front of or behind the receiving 
slit as illustrated in Fig. 14 for the asymmetric diffiac- 
tion case. 

The aberration profile Jdf(^) for this condition is a 
impulse function of angular width A^^ = DIR (rad) 
where D is the width of the defocussed beam at the 
receiving slit. Assuming the equatorial divergence a is 
small (ie., 2° or smaller) the angular width zl^p of the 
aberration profile is given by 



D 



R^-R 



R 



(13) 



Under asymmetric diffraction conditions R2 and R are 
related by the equation. 



R^=R 



sin(2ft-6)) 
sin6> 



(14) 



SO that the width A^p of the aberration function is. 



A^^=a 



sin(2ft-6>) 



sm6> 



(15) 



In asymmetric diffiaction, the breadth of diffraction 
lines increases as the deviation from the symmetric 
condition, \(0-Q\, increases. Defocussing is larger at 
low 2ft angles and varies more rapidly with 6>- ft at 
low 2ft values. Conversely at high 2ft values, diffiac- 
tometers will tolerate quite large errors in 6>- ft(ie., up 
to ±5° at 2ft- 150°) without the effects of defocussing 
being detectable in the line breadth. Also, by reducing 
the angle of divergence a the defocussing A^^ can also 
be reduced, but at the expense of the diffiacted intensi- 
ty. The plot of zIdf vs {(O- ft) for a range of 2ft angles 
from 30° to 150° is shown in Fig. 15. 




29.90 



29.95 20 30. 00 




-[1 = 50 cm^ 
■ ^i = 100 cm ' 
- ^i = 200 cm ■ 



89.85 



89.90 



89.95 20 90.00 



Fig. 13. Specimen transparency aberration profiles at 20= 30° and 90° for linear attenu- 
ation coefficients of 50, 100, and 200 cm~\ 
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JDF(e) 



M ► 

Adf = D/R 



X-ray Source 



Fig. 14. Defocussing under asymmetric diffraction conditions giving an impulse shaped aberra- 
tion function /opCe) of angular width ^^jp. 



20 20 

150° 120° 90° 30° 30° 60° 90° 120° 150° 




Fig. 15. The variation of defocussing width ^^p with (co- 9)^t different 26 
diffraction angles. 



4.6 Axial Divergence 

In a laboratory diffractometer only a small fraction 
of the photons that form the incident beam emerge from 
an x-ray source parallel to the equatorial plane. 
Likewise most of the x rays that reach the detector slit 
after diffracting from the sample are angled to the equa- 
torial plane. Under these circumstances, a diffractome- 
ter will record x-ray counts over a range of angles 2<^ 
other than the true diffraction angle 20. The only rays 
for which 2<^ = 2^ will be those propagating parallel to 
the equatorial plane and incident on the diffractometer 



axis. In practice, axial divergence is most readily recog- 
nised by the asymmetry it introduces into low angle dif- 
fraction lines (2^<30'^) where the low angle tails 
extend further than the high angles tails. 

For a particular ray path the measured diffraction 
angle 2 <^ for a true diffraction angle 20 depends on the 
axial divergence p and y in the incident and diffracted 
rays (see Fig. 16). Assuming small angles for p and % 
the difference £=2<^-2^is given by, 

£ = 2(l)-20 = l5rcosQc20 -^— ^cot2a (16) 
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Incident 
Ray 



Fig. 16. An illustration of axial divergence in which the incident ray and diffracted ray are at 
angles of j3 and /relative to the equatorial plane. Although the diffraction angle at the sample is 
29, the diffracted beam is recorded at an angle 2<^. 



From this equation it is evident that the effect of axial 
divergence on e will not only be large when 2^ is small, 
but also when 20 is large (ie., 20> 150°) where the 
asymmetry is opposite to that at low 20. Although axial 
divergence is as strong in very high angle lines as it is 
in low angle lines, the positive asymmetry developed 
tends to be less noticeable as it is overshadowed by the 
dispersion of the emission profile. 

When axial divergence in the incident beam is small 
(P ~ 0) as in many synchrotron systems or laboratory 
diffractometers with very narrow incident beam Soller 
slits, all the axial divergence arises from the diffracted 
beam and Eq. (16) reduces to. 



£=—^cot20, 

2 



(16a) 



For this condition, axial divergence is absent at 20 = 
90°, £ < when 20 < 90° and £ > when 20 > 90°. The 
aberration function Jax(£) for this condition can be 
derived analytically [22,23,3,4] by considering an axi- 
ally parallel beam from a line source incident on a nar- 
row capillary specimen of randomly oriented crystal- 
lites as shown below in Fig. 1 7. 




specimen 



Fig. 17. Arcs formed in the vicinity of the receiving slit by the dif- 
fraction cones radiating from all points across the sample axis. The 
axial divergence aberration function is the profile shape detected by 
an infinitely narrow receiving slit as it scans through the arcs. 

The emergent diffracted beam from each individual 
ray in the incident beam in Fig. 18 is a collection of 
radiating cones with a semi-angle 20 and Jax(^) can be 
derived by determining the intensity profile when the 
receiving slit is scanned across these cones. The param- 
eters of the diffractometer that define Jax(^) in this 
instance are Lg, the axial length of the specimen bathed 
in X rays and the length L, of the receiving slit and the 
function Jax(^) is given by, 



^Ax(^)-- 



^Ax(^) 



is is 



for£={0,f;} (17a) 



for e ={s^, e^} (17b) 



where e, = -_ — — \ rad 



and 



£, =- 



2 [ 2R 
cot20 fL^+L^ 



2 I 27? 



rad. 



Examples of this aberration function at 20= 10° and 
50° are shown in Fig. 1 8 for conditions that resemble a 
laboratory diffractometer with very narrow incident 
beam Soller slits (ie., <1°), R = 215 mm, an illuminat- 
ed specimen length Ls = 12 mm and a receiving slit of 
length Lr = 1 6 mm. 

When allowance is made for axial divergence in both 
the incident and diffracted beams and for the presence 
of Soller slits in each of these beams, the calculation of 
the aberration function Jax(^) can no longer be done 
analytically. Eastabrook [24] first demonstrated how to 
calculate the axial divergence aberration profile for a 
conventional diffractometer, but restricted his discus- 
sion to instrumental conditions that could be solved 
analytically rather than for conditions that were widely 
used in practice. Pike [25] generalised the application 
conditions and included the effect of Soller slits, but 
restricted his calculations to the determination of the 
centre of gravity and variance of profiles rather than the 
aberration profile itself More recently Cheary and 
Coelho [3,4] have developed a semi-analytical 
approach to the calculation and their results have been 
incorporated into a profile refinement procedure. In 
short, their procedure consists of, 
(i) calculating the analytical aberration function 
Jax(A^) arising from incident rays all with the 
same axial divergence p. The instrument parame- 
ters required to define this function are the axial 
lengths L^, Lg, and L, of the x-ray source, the sam- 
ple and the receiving slit, 
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Fig. 18. Axial divergence aberration profiles Jax(^) for a diffractometer with no axial divergence 
in the incident beam. This was calculated using Eq. (17) for Lg = 10 mm, Li^= 15 mm, and 7? = 
215 mm at 20= 10° and 29= 50°. A very narrow Lorentzian profile was convoluted with J^^(£) 
to overcome the infinities at £= 0. 



(ii) incorporating the Soller slits into the calculation as 
angular intensity filters on the axial divergence P 
and yin the incident and diffracted beams, respec- 
tively The transmission functions, S^(P) and Sj)(f), 
for the incident and diffracted beam Soller slits, 
respectively are each triangle functions with 
100 % transmission at ^8= and 7= 0, and % 
transmission at j8 = ± A^/l and 7= ± A^^/l where A^ 
and z\d are the angular apertures of the slits (see 
Fig. 9 for definition of z\), 

(iii) calculate the full aberration profile /axC^) by inte- 
grating the aberration functions determined at 
each P across all the allowed P values, 

J^ie)= j J^iP,e)S,ip)S^i^p,e\)dp (18) 

In practice this integration is converted to a sum- 
mation, but great care is needed carrying out the 
summation because of singularities in the AxCA^) 
functions. Considerable savings in computing 
time without loss of accuracy can be achieved by 
using a smoothing procedure in which the function 
Jax(£) is convoluted in e space with an impulse of 
width = step size of the data being fitted. The sav- 
ing in computing time comes from a reduction in 
the number of summation points at discrete values 
of P required to accurately reproduce Eq. (18) 
using a summation. 
In the absence of Soller slits the maximum axial 
divergence in the incident and diffracted beams, A'^ 
and A^'''' respectively is determined by the length L^, 
Zg, and L, and given by. 



tanzi^ 



R 



and tan zip 



R 



(19) 



Under these circumstances the axial divergence can be 
as large as 10° in the incident beam and 10"^ in the dif- 
fracted beam for a diffractometer with a 12 mm long x- 
ray source, a 12 mm receiving slit and a 20 mm wide 
specimen. However, most commercial diffractometers 
are supplied with Soller slits having angular apertures 
A between 2° and 5°, and in such cases the breadth of 
the aberration profile is reduced quite dramatically. 
This is illustrated in Fig. 19 which shows aberration 
profile shapes calculated using the method of Cheary 
and Coelho for diffractometer configurations that 
include Soller slits as well as a configuration with no 
Soller slits. When narrow Soller slits are included in the 
beam path (ie., A^ and A^ ^ 2°), the dimensions L^, L,, 
and L^ tend to be redundant parameters in the calcula- 
tion of Jax(£) as the Soller slits control the maximum 
axial divergence in the incident and diffracted beams. 

The major change brought on by axial divergence is 
extension of the low angle tails of profiles below 2^~ 
50"^. At very high angles, 20> 150°, the asymmetry 
reverses and the extension of the high angle tails 
increases with increasing 20. Broadening from axial 
divergence is evident at all 20 values but passes 
through a minimum in the region 20 ~ 110°. However, 
the shift in peak angle 20,^^ at I,^^ relative to the true 
2^ is small at all 2 lvalues unlike the shift in centre of 
gravity 20^^ -20q which varies quite considerably from 
very large negative values, near 2^= 0°, to very large 
positive values as 2^ approaches 180° [12]. 
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Fig. 19. Axial divergence aberration profiles /ax(^) calculated for various combinations of inci- 
dent beam and diffracted beam Soller slits (or no Soller slits) 29= 20° and 29= 90°. Instrument 
conditions: Diffractometer radius 215 mm, L^= 12 mm, L^ = 25 mm, L^= 12 mm. The three 
combination used are (i) No Soller slits, (ii) 2° Soller slits in the incident beam only and, (iii) 2° 
Soller slits in both the incident and diffracted beams. All the aberration profiles are normalised to 
the same L^^^. 



In practice it is not always possible to calculate the 
exact form of the axial divergence aberration function 
for a particular specimen/diffractometer configuration. 
The two main reasons for this are, 
(i) in specimens with strong preferred orientation, 
such as thin films and rolled or extruded metals, 
the diffraction cones are no longer of uniform 
intensity along the arcs of the diffraction cones, 
(ii) the inclusion of a monochromator in the beam 
path reduces the axial divergence. When either 
diffracted beam or an incident beam monochroma- 
tor is present the optical path length of the beam is 
extended. With a diffracted beam monochromator 
the optical path length of the diffracted beam is 
extended by the optical path between the receiving 
slit and the detector slit. In a diffractometer with a 
graphite monochromator tuned to Cu Ka radiation 
the optical path length in the diffracted beam is 
increased by ^100 mm giving an effective radius 
of the detector arm of 100 + 7? mm. An example of 
the extent to which axial divergence is reduced by 
introducing a graphite diffracted beam monochro- 
mator is illustrated in Fig. 20. Incident beam 
monochromators also increase the optical path 
length of the incident beam. For the most common 
type of monochromator in use, the asymmetrical- 
ly Ge ground and bent monochromator, the path is 
increased by over 300 mm. Monochromators also 
act as angular intensity filters and their effect on a 
profile is similar to the addition of Soller slits in 
the beam path [26]. The effect of a monochroma- 
tor can therefore be represented as a Soller slit in 
a profile fitting model. 



no 
monochromator 




with graphite 
monochromator 



20.9 



— I — 
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Fig. 20. Effect on the 100 profile from LaB^ powder of including a 
graphite monochromator in the diffracted beam in a diffractometer 
with no diffracted beam Soller slits. The two profiles were nor- 
malised to the same peak intensity Z^^^- 



4.7 Linear Position Sensitive Detector (LPSD) 
Aberrations 

LPS detectors with angular windows up to 10° (at 
200 mm) are used in reflection mode on commercial 
diffractometers to increase the data collection rate par- 
ticularly for kinetic based studies. In most cases LPS 
detectors are used in stationary mode and only a fixed 
angular region of a pattern is recorded. A number of 
manufacturers offer LPSD systems where the detector 
can be scanned in 0-20 mode with its centre main- 
tained at the focus condition. As the detector is scanned 
the total diffraction pattern is formed by adding togeth- 
er and averaging the patterns recorded at each step. In 
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this mode a diffraction pattern can be accumulated very 
rapidly with excellent counting statistics because indi- 
vidual lines are within the detector window for a 
considerable time. For example a LPSD with an angu- 
lar window of 10° moving at 5° per min corresponds 
to a diffraction peak being detected for 120 s. 
Consequently, in this mode, a 100° 20 diffraction pat- 
tern may only take approximately 20 min to collect. To 
obtain the same level of counting precision in a conven- 
tional single slit diffractometer with 0.02° 20 steps 
would take approximately one week. In practice, this 
gain in time is never fully realised because of the longer 
dead times of LPSDs and the damage that can be 
caused on the anode wire by very high localised count 
rates. Also the diffraction peaks recorded at off-centre 
positions along the detector window are broadened and 
asymmetric. Some manufacturers make the angular 
window a software adjustable parameter to maintain 
the resolution although this means that only a fraction 
of the detector is being used. 

In LPSD systems, the receiving slit aberration is no 
longer relevant and the flat specimen aberration is 
replaced by an aberration function that embodies three 
effects which are folded together in the final function; 

• flat specimen error including defocussing, 

• parallax error, 

• thermal noise. 

A full treatment of these aberrations and how they are 
modified by a scanning LPSD system is given in 
Cheary and Coelho [14]. For a stationary LPSD the 
aberrations are discussed in Sees. 4.7.1-4.7.3. 

4.7.1 LPSD Flat Specimen Error /psd(£) Including 
Defocussing 

In a conventional diffractometer defocussing and flat 
specimen error can be convoluted together independ- 



ently, but this is no longer valid in a LPSD system. The 
effect of both the specimen and the detector being flat 
is that the defocussing is no longer symmetric about the 
centre of the LPSD. In a stationary LPSD the aberration 
function Jpsd(^) depends on the 20 value of a peak and 
on the offset angle P of the recorded peak from the cen- 
tre of the LPSD. The angular variables used to define 
the LPSD system are given in Fig. 21 where the ^ is 
limited by the angle of divergence a of the incident 
beam (ie., -a/2 < f < +a/2). 

When the diffractometer is operated symmetrically 
and the incident beam is centred on the diffractometer 
axis then assuming small angles for e and f, the differ- 
ence 8=2(1) -20 is related parabolically to f, 



e = aC-bC or iC+Q' 



where a = cos p 



s'miO + p/2 
sm(0-p/2 



-1, b-- 



("20) (20) 



71 sin2^cos P 
\S0sm\0-p/2) 



> 0, £f. = — and Cn = — - The equation for the aberra- 

'4b '2b 

tion function Jpsd(^) is obtained by transforming the 
intensity across the incident beam IB(^ from f space 
into £ space, ie., 



Jp3Je) = IB(C) 



dC 



de 



(21) 



As the incident beam intensity IB(^ is reasonably con- 
stant then the un-normahsed form of Jpsd(£) = |d^/d£ |. 
This is readily calculated by differentiating the trans- 
formation Eq. (20) and has the form, 



^PSD \^) ' 



K 



V^T-^ 



(22) 



LPSD 




Fig, 21, Beam geometry for an LPSD diffractometer showing the various angular 
parameters used in the derivation of the aberration function JpsoC^)- 
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where £<£q. When calculating this function allowance 
has to be made for the fact that Eq. (20) can be double 
valued. Physically this means that more than one part 
of the incident beam contributes at a particular £ value 
and Jpsd(^) possesses a discontinuity at the boundary 
between two rays and one ray contributing to the aber- 
ration function. This is illustrated in Fig. 22 for profile 
arising from a central ray entering the LPSD at 0.5° off- 
centre when set at a low 20 (ie., j3 = 0.5°, 2^=20°). 
When £q is outside the range of £ values dictated by the 
limiting values of ^ in the incident beam (typically 
±0.5°) then there are no discontinuities within Jpsd(^)? 
no infinities and Jpsd(^) is finite at all £ as illustrated in 
Fig. 23 for the profile with (5 = -2°. It should be noted 
that the aberration functions for a particular 2^ are not 
symmetric with respect to (5. In Fig. 22 the profiles for 
P = +2° and P = -2° at20= 20° are distinctly different 
and have different limits in £. 



p=+2' 




riQ) 



Fig. 22. LPSD flat specimen error/defocussing aberration function 
Jpsd(£) for an incident beam of divergence a=l° at 20 = 20°. 
Examples are given of Jpsd(£) near the centre of tlie PSD (j8= 0.5) 
and at 2° either side of centre (ie., P = +2° and -2°). 
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X-rays Incident on Detector 

Fig. 23. Aberration profiles for parallax error and the effect of the 
absorption jUg of the detector gas. 



4.7.2 Parallax Error 

When a diffracted x-ray photon enters a LPSD its 
path in the detector gas is not perpendicular to the 
anode wire, except at the centre of the detector, and 
additional broadening, known as parallax broadening, 
is introduced into off-centre diffraction lines. This aris- 
es because the ionisation caused by incoming photons 
is likely to occur at any point along its path, but the sub- 
sequent avalanche is perpendicular to the anode wire 
(see Fig. 24). For a detector gas with low absorption, 
the profile shape recorded by the LPSD at an angle ±(5 
from the centre of the LPSD is a unit area impulse func- 
tion with a width zipx (°20) given by 



JP\D 
R 



(23) 



where D is the depth of the detector and R is the radius 
of the diffractometer. Parallax error can be reduced 
considerably by increasing the detector gas pressure. In 
this way the x rays are absorbed by a thin layer of gas 
beneath the detector window, in which case the diffrac- 
tometer radius R is defined by the front window of the 
detector rather than its anode wire. 




Fig. 24. Ka^Oj emission profile fiom copper recorded using the 400 
line from a sihcon single crystal. 

Most detectors are designed with a quantum efficien- 
cy of at least 80 % and have depths somewhere 
between 5 mm and 1 mm. The aberration profiles for 
these conditions resemble those shown in Fig. 23 where 
the absorption of the detector fi^^^ > 0. The shape of the 
aberration profiles for parallax broadening J^psd(^) at 
an angle P to the centre of the LPSD is given by. 



^uPSD \^) ~ 



/I,/? 



exp — 






|)8|[l-exp(-/iO)] 



(24) 
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where jU^ is the hnear attenuation coefficient of the 
detector gas [27]. Parallax broadening can be quite 
large for long detectors and is usually the dominant 
aberration when LPSDs are used for peaks above 
2^-40°. For example, when ^8 = ±5° and the detector 
has a depth Z) ~ 8 mm, the breadth of the parallax func- 
tion Ap can be as large as 0.2°. 

4.7.3 Thermal Noise 

The ultimate resolution of a LPSD is the spatial 
uncertainty of the position measurement of an individ- 
ual x-ray photon incident normal to the anode wire of 
the detector. This is controlled by the spatial broaden- 
ing of the discharge produced by incoming photons, by 
the thermal noise generated in the anode wire and by 
the accuracy of the position-encoding electronics [28]. 
In most instances the thermal noise profile Jnpsd(^) is 
modelled as a Gaussian function with a FWHM a, 

TN(£) = (2/cr)(ln(2;r)'^'exp[-ln(2)(2 e/ a)2] (25) 

This function represents the distribution of "measured 
positions" of a discharge caused by an incoming x-ray 
photon when the photon is incident normally at exactly 
the same position along the wire. In most LPSDs, the 
angular width a is equivalent to a positional uncertain- 
ty along the anode wire between Ax ^ 0.04 mm and 
0.20 mm where a=Ax/R radians. This effect con- 
tributes to the broadening much like that from the typ- 
ical receiving slit width in a conventional diffractome- 
ter. Smaller Ax values down to 0.040 mm are possible 
when the gas in the LPSD is under pressure. Most 
LPSDs have a positional resolution Ax around 0.1 mm 
although those that utilise delay lines can be as large 
0.2 mm. 



5. The Wavelength Distribution in 
Laboratory Diffractometers 

The natural shape of the energy distribution W(E) of 
a single characteristic x-ray emission line well above 
the threshold energy is Lorentzian and given by 



WiE) = 



nin 



{E-E,f+{ri2f 



(26) 



where £"0 is the peak emission energy and Tis the life- 
time broadening given by the sum of the widths of the 
two relevant atomic levels involved in the transition. 



When a Lorentzian energy spectrum is transformed into 
20 space it remains Lorentzian, provided dE'/d2^ is rel- 
atively constant across the profile, with a FWHM F20 in 
2^ space, given by 



2rtanei80 



n 



(°20) 



(27) 



where Fand £"0 are in eV. The form of the Kaia2 emis- 
sion profile from Cu is shown in Fig. 24. For all of the 
transition element anodes used in x-ray diffraction nei- 
ther the Kai line nor the Ka2 line has a Lorentzian 
shape and both the Kai and Ka2 lines are asymmetric 
with extended high angle tails. Moreover, the asymme- 
tries and FWHM values of the Kai and Ka2 peaks are 
different. This is particularly evident with the first 
series of transition element and is related to the anom- 
aly in the atomic number dependence of the atomic 
level widths of the Lh and Lm levels obtained from x- 
ray photoelectron spectroscopy data and x-ray emission 
spectroscopy [29]. The line widths F, asymmetry 
indices /cand energies E^ of the four transition elements 
commonly used for x-ray targets are given in Table 1 . 

Table 1. Line widths F, asymmetry indices ?cand emission energies 
Eq for tlie Ka^ and Ka2 emissions of selected transition elements 
(data from Salem and Lee [30]) 







Kai 






Ka2 




Element 


r(eV) 


K 


£o(keV) 


r(eV) 


K 


£o(keV) 


Cr 


2.16 


1.38 


5.415 


2.75 


1.18 


5.406 


Fe 


2.35 


1.43 


6.404 


2.84 


1.25 


6.391 


Co 


2.87 


1.32 


6.930 


3.59 


1.25 


6.915 


Cu 


2.56 


1.12 


8.048 


4.05 


1.10 


8.028 



In FPPF an accurate model of the emission profile is 
essential particularly for the analysis of high angle pro- 
files. The FWHM r20of the Kai line from Cu and other 
commonly used transition metal anodes is shown in 
Fig. 25. At 2e = 40° the FWHM F^e is less than or equal 
to 0.01° {20) is considerably smaller than the FWHM 
values of actual diffraction lines from commercial dif- 
fractometers which are typically in the range 0.07° to 
0.10° 20 [21] depending on the choice of slits. In this 
region the contribution of the emission profile is 
swamped by the geometrical aberrations. When 20 > 
60°, however, the emission profile tends to dominate 
over the geometrical aberrations. Once 20> 100°, the 
total breadth of the geometrical aberrations is relative- 
ly minor and the profile shape conforms closely to the 
emission profile. 



Volume 1 09, Number 1 , January-February 2004 

Journal of Research of the National Institute of Standards and Technology 



0.15 




Fig. 25. Full width at half maximum intensity of the Ka^ Hne i~20as 
a function of 29 for four different target materials. 



The natural asymmetry of the emission hnes arises 
from the multiplet structure of the transitions. In addi- 
tion to the main transitions involving the change in 
vacancy state 1^ -^ 2p, it has been recognised that ?>d 
spectator transitions also contribute up to 30 % of the 
Kaia2 emissions [31]. In these transitions the atom is 
doubly ionised and the actual vacancy transition is still 
\s ^2p, but the second vacancy in the 3d level is not 
directly involved in the transition. The notation for this 
transition is \sM -^ 2p3d. A phenomenological repre- 
sentation for accurately describing the asymmetric Cu 
emission profile was first used by Berger [32]. In this 
model the Kai and Ka2 lines were each represented as 
two Lorentzian profiles as shown in Fig. 26. This rep- 
resentation was also used successfiilly on accurate 
spectroscopic data from Cu obtained by Hartwig et al. 
[33]. A systematic study of the Ka and K(3 emission 
profiles from Cr, Mn, Fe, Co, Ni, and Cu [34] has 
shown that the phenomenological representation can be 




1.538 1.540 1.542 1.544 1.546 

Wavelength (A) 

Fig. 26. Phenomenological representation of the Cu Ka emission 
profile based on four Lorentzians. 



used to accurately represent these elements down to an 
R factor of 1 %, although in some cases it is necessary 
to use up to seven Lorentzians. 

Another feature of Ka emission lines which needs to 
be included in an accurate line profile fitting model is 
the satellite multiplet structure in the high energy tail as 
shown in Fig. 27. As a group these have an intensity of 
-0.6 % of the Kai emission line in the case of Cu ris- 
ing uniformly with decreasing atomic number up to 
-1.4 % for Cr [35]. Satellites lines are also evident in 
KPiPg spectra and appear on both the low energy and 
high energy tails as the K(3' and K(3" lines [36]. The 
Ka3a4 non-diagram lines arise from the transition 
ls2p—>2p^ in which the actual vacancy transition is 
Is ^2p, but in the presence of a 2p spectator hole 
[37,38]. For most of the first transition series the Ka 
satellite structure can be fitted accurately with four or 
five Lorentzians [35,38]. In most x-ray diffraction stud- 
ies this level of precision in fitting the satellites is 
unnecessary. For Cu it is sufficient to represent the Ka 
satellite group as a single broad Lorentzian so that the 
total Ka spectrum can be represented by five 
Lorentzians as given in Table 2. 



FWHM Ka. -0.045^26 




Fig. 27. Cu Ka emission profile showing the sateUite group of lines 
and the extent of the tails from the Kai and Ka2 emission lines. This 
profile was recorded using the 400 line from a silicon single crystal 
wafer. 

Table 2. Relative intensities /q (areas), wavelengths A and lifetimes 
widths Fx (in A units) for representing the Cu Ka spectrum by five 
Lorentzians. The Kaia2 data were taken from Holtzer et al. [34] and 
the satellite data from Cheary and Coelho [2]. 



Emission 


A(A) 


Relative /q 


r^xio^A 


line 








Ka,3 


1.540591 


0.5710 


0.437 


Ka,b 


1.541064 


0.0789 


0.643 


Kaja 


1.544399 


0.2328 


0.513 


Kajb 


1.544686 


0.1036 


0.687 


Ka3a4 Satellites 


1.534753 


0.0137 


3.686 
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In the transformation of the distribution W{X) from 
wavelength space to 20 space at 2^< 130"^ the 
Lorentzian shape in A maps directly into a Lorentzian 
in 2^ with minimal error. At higher angles the effects of 
dispersion between X and 2^ become increasingly evi- 
dent [12] and the transformation needs to be carried 
more accurately The frill expression for the trans- 
formed wavelength distribution in 2^ space for a sum 
of Lorentzians in A space is 



WilO) 



1 



dA 



dcoseti[{rj,,/2f+{2dsme-Xf] 



(28) 



where /q/ and ^;^J are the relative areas and wavelength 
FWHM of the /th Lorentzian. The d spacing is defined 
in relation to the reference wavelength X^f of the emis- 
sion line which in the TOPAS implementation of this 
procedure is the wavelength of the highest intensity 
Lorentzian. In the case of the Cu Ka spectrum given in 
Table 2, the reference wavelength is that of the Kai^ 
component (A =1.540591 A) so that X,Qf= X^^i^ = 
IdsinO^^^^, At very high 20 values (ie., 20 > 150°) the 
cosO terms arising from d20/dA can change quite sig- 
nificantly over a profile and elevate the intensity in the 
high angle tail. In addition to dispersion, the distortion- 
al effects of the Lorentz factor on high profiles should 
also be incorporated when analysing profile shapes at 
very high angles [39,40,41]. 



\if 







As the natural shape of an emission profile is 
Lorentzian, the tails can extend a considerable distance 
from the central peak as shown in Fig. 27. Most diffrac- 
tometers, however, operate with K(3 filtering or some 
form of monochromatisation and as such the tails are 
attenuated to varying degrees. When a Ni K(3 is includ- 
ed in the beamline of Cu Ka instrument, the attenuation 
appears to be more or less uniform across the profile 
except below the Ni K absorption edge. Careful analy- 
sis however, reveals a small variation in the attenuation 
across the profile owing to the increase in linear atten- 
uation coefficient with increasing wavelength (ie., 
|i a ^^), but this is generally a small effect and is not 
expected to affect the profile shape significantly [12]. 
The inclusion of a curved graphite monochromator in 
the diffracted beam of a diffractometer greatly reduces 
the range of wavelengths entering the detector resulting 
in profile tails that diminish more rapidly than the nat- 
ural emission profile as shown in Fig. 28. These mono- 
chromators can also affect the relative intensity /Ka2/^Kai 
ratio by up to ±10 % depending on the alignment and 
setting of the crystal. For example, in the monochro- 
mated spectrum shown in Fig. 28 the relative intensity 
of the CuKai/Ka2 is approximately 0.46 rather than 
0.50 as in the unfiltered spectrum. 

Commercial pyrolytic graphite monochromators are 
strongly oriented polycrystals bent by a hot pressing 
operation. The focussing is imperfect and the resolution 
is relatively poor. Better resolution and, in some cases, 
high intensities can be obtained from "ground and 



Cu Ka Unfiltered 




KffccE of 

Graphite 

MurtoL-hrumaior 



70 



Fig. 28. Cu Ka emission profile obtained using the 400 line from a silicon single crystal 
wafer. Each pattern was recorded sequentially using the same sample, first with no filter or 
monochromator in the beamline, then with a NiK(3 filter, and finally with a standard curved 
graphite diffracted beam monochromator. 
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bent" (Johannson) single crystals. These crystals 
achieve perfect focussing when correctly aligned and 
are able to select a very narrow wavelength band. The 
most common materials used for Johannson monochro- 
mators on conventional diffractometers are quartz, ger- 
manium and silicon. The example shown in Fig. 29 
shows the wavelength spectrum from an asymmetrical- 
ly cut ground and bent Ge crystal used as an incident 
beam monochromator. The wavelength passband is 
narrow enough to remove 99.98 % of the Ka2 compo- 
nent and 1 00 % Ka satellites from the CuKa spectrum 
and, almost completely eradicate the Lorentzian tails of 
the emission profile. 

In the presence of monochromators, even low reso- 
lution graphite monochromators, it is no longer possi- 
ble to accurately model the wavelength distribution 
from an x-ray tube using the tabulated unfiltered spec- 
tra such as the one in Table 2. Although first principles 
calculations of the wavelength transmission fiinction 
through ideal monochromator are possible, it is current- 
ly more practical to determine experimentally a 
"learned" spectrum for the curved graphite, germani- 
um, quartz and lithium fluoride monochromators used 
in the majority of monochromated laboratory powder 
diffractometers. This can be done by modifying the 
"sum of Lorentzians" representation in energy space or 
X space to fit the spectrum entering the detector. In 
broad terms, monochromators reduce the width of the 
wavelength distribution and tend to truncate the tails of 
spectra. A number of approaches can be used to accom- 
modate these changes, 

• represent the components of the wavelength distribu- 
tion as Voigt or pseudo-Voigt fiinctions rather than 



Lorentzians, to limit the extension of the profile tails, 
and modify the relative intensities of the component. 
• represent the effect of the monochromator as a wave- 
length filter with a transmission function T{X) repre- 
sented by a simple function, such as a split pseudo- 
Voigt or split Pearson VII function with up to four 
refineable parameters, which operates on the tabulat- 
ed Lorentzian emission profile data. The split func- 
tions are used to incorporate asymmetry in the T{X). 
In both cases the parameters of each representation are 
obtained by analysing and fitting the high angle profiles 
from either a reference line profile standard, such as 
LaBg SRM 660a, or single crystal disc such as a 111 
wafer of silicon. 

The inclusion of a parabolic multilayer mirror in the 
incident beam of a diffractometer can also introduce a 
distortion into the wavelength spectrum [42] . This hap- 
pens because the Kai and Ka2 components of the spec- 
trum reflect off the mirror in slightly different direc- 
tions as shown in Fig. 31. The separation of the Kai 
and Ka2 peak maxima, A2^Ka2i = 2^ (Ka2) - 2^ (KaJ, 
in a line profile is then either larger or smaller than the 
same profile from a diffractometer with no mirror. 
When a diffractometer is set up with its mirror in the 
orientation shown in Fig. 30, the separation A2^Ktx2i is 
smaller than the value expected from the known Kai 
and Ka2 wavelengths by an amount corresponding to 
the difference Ai//^in the directions of the two incident 
beams on to the specimen. According to Toraya and 
Hibina [42] the difference 8^ for a conventional long 
fine focus x-ray with a projected target width of 0.04 
mm can rise to 0.0017° for a high resolution mirror, but 
decreases to a negligible level for low resolution mir- 
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Fig. 29. Wavelength spectrum emerging from an asymmetrically cut Gel 1 1 ground and bent incident beam monochromator presented lin- 
early and logarithmically. The Ka satellites are completely removed, but the Ysl^ is still present at -0.02 % of Ka^ even in a well aligned 
system. 
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Fig. 30. Reflections of Ka^ and Ka2 wavelengths from a parabolic multilayer 
mirror diffracting of a powder specimen. 
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Fig. 31. Results of fitting HI and 220 profiles from MgO with the parameters shown in Table 3. 



rors. In any high accuracy lattice parameter determina- 
tion using this type of diffractometer it is necessary to 
incorporate the wavelength dependent zero error into 
the analysis. 



6. Fundamental Parameters Profile 
Fitting (FPPF) in Practice 

In practice FPPF requires accurate numerical proce- 
dures for carrying out multiple convolution integrals. 
This can be done by representing the calculated profiles 
as a histogram and reducing the convolution integral to 
a summation. To avoid systematic errors with this 
approach the angular step size between calculated 
intensities needs to be very small but in so doing the 
operation becomes very time consuming [2]. In the 
XFIT and TOPAS implementations of FPPF a semi- 
analytical procedure has been developed for convolut- 
ing the aberration functions. In this procedure the two 
aberration functions being folded together are calculat- 
ed at the same 20 values as the measured data points 
and then a continuous function is formed by interpolat- 
ing between the calculated points. As the calculated 
functions are then each a series of linear sections it is 
possible to calculate the convolution integral analyti- 
cally [2]. Some difficulties are experienced with aber- 
ration functions that possess singularities, but in all 



cases the functions are integrable and the effect of the 
singularity can be overcome either by a convolution 
process or by a smoothing operation [4]. 

Although FPPF is a powerful method of profile 
analysis it can fail when used without recognising its 
limitations. The physical parameters describing a dif- 
fractometer are not statistically independent within the 
least squares refinement procedure and strong correla- 
tions exist between many of the refineable parameters. 
For example, when the angular acceptance angles, zij 
and z\d? of the incident beam and diffracted beam Soller 
slits are refined independently to a set of profiles, vari- 
ous combinations of refined values zij and zip can be 
obtained that give equally good fits to the data. The rea- 
son for this is that these two parameters are very strong- 
ly correlated and the values can in fact be interchanged 
without changing either the shape of a synthesised pro- 
file significantly or the quality of the fit. As a conse- 
quence the best fit with an unconstrained refinement 
often corresponds to refined parameters that do not 
make physical sense. When first using the FPPF 
approach it is advisable to investigate the validity of the 
refined parameters for a particular diffractometer using 
a well crystallised reference specimen and establish the 
parameters of the instrument that allow profiles over 
the whole 2 Grange to be defined with one set of values. 
It would be unrealistic to expect the refined parameters 
for a diffractometer to match the directly measured val- 
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ues exactly as there are a number of second order 
effects in diffractometer profiles that are not incorpo- 
rated in the fitting model. Moreover, not all of the 
instrumental aberrations are independent and as such 
the convolution model is not strictly valid with certain 
combinations of aberrations [8]. Nevertheless, experi- 
ence with most of the commercially available diffrac- 
tometers has shown that refined values reasonably 
close to the true instrumental values can be expected. In 
the original investigation of fundamental parameters 
fitting [2] it was shown that when deliberate changes 
were made to the diffractometer set-up, such as chang- 
ing the receiving slit length, receiving slit width and 
diffractometer radii, the change in the refined instru- 
ment values corresponded well with the actual changes 
despite the limitations in the axial divergence model 
used at that time. In general, therefore, if physically 
unrealistic instrumental parameters are required to 
describe the diffraction pattern of a reference material 
then there is either a deficiency in the model used to 
describe the diffractometer, some sort of mis-setting of 
the diffractometer or the refinement is trapped in a false 
minimum. 

Once the refined instrument parameters have been 
established (ie., the dimensions and apertures of the 
various slits, the source size and the wavelength distri- 
bution), and these are in reasonable agreement with the 
actual values, the only instrument parameter that may 
need refinement from specimen to specimen is the lin- 
ear attenuation coefficient //. Absorption effects in a 
line profile are only evident when // < 1 00 cm"^ and 
under these circumstances there is some justification in 
making // a refineable parameter. In specimens with 
// > 200 cm"^ refinement of ji normally has little effect 
on the profile shape because the specimen transparency 
profile is relatively narrow at all 20 angles (ie., 
<0.01°26*). Consequently, // is fixed at a representative 
value and not refined. 

Even in well crystallised powders with crystallite 
sizes up to 2 |xm crystallite size broadening is 
detectable in high angle lines. Conversely, in many 
powders it is not uncommon to have crystallite sizes 
down to 0.5 |xm and crystallite size broadening is evi- 
dent even in low angle lines. Also, in many crushed 
powders the action of crushing even in a standard pes- 
tle and mortar can induce microstrain broadening. In 
TOPAS and XFIT therefore, the apparent crystallite 
size r^pp and the percent microstrain £^ are frequently 
included as refineable parameters although their contri- 
bution is generally small. Various profile functions can 
be adopted for crystallite size and microstrain and com- 
monly take the form of a Lorentzian function or a 



Gaussian function. The most common unit area func- 
tions used for crystallite size and microstrain, B^^^(20) 
and B^(20), respectively, are 



K.A20)-. 



H^J2k 



{2e-2e,fHH^j2f 



BA20)- 



H.. 



exp 



-In 2 



^2(20-20^^ 



(29a) 



(29b) 



where the FWHM of the Lorentzian = H^^ = 180 
yU(7iT,^ cosQ, ) °2e md the FWHM of the Gaussian = 
H^ = (18£uV21n2/5;r) tan ^o °2a In some crushed 
powders it has been necessary to use a Lorentzian func- 
tion to represent the microstrain but with a FWHM = 
TjimiO where r/ is a constant related to the microstrain 

[4]. 

The quality of the fits obtainable from the FPPF 
approach and the degree to which the refined instru- 
ment parameters agree with the actual values is illus- 
trated below for data collected from a polycrystalline 
MgO reference specimen over the range 20 = 36° to 
150° using CuKa radiation. The specimen used in this 
instance was a 20 mm disc prepared by sintering at 
1400 °C for 24 h in air. When a small segment of the 
disc was examined in a scanning electron microscope 
the vast majority of crystallites were between 1 |xm and 
2 |xm in diameter. Table 3 shows the refined values 
obtained for the MgO data from TOPAS whilst Fig. 30 
shows the fits obtained to two of the profiles that are 
sensitive to the instrument parameters and the specimen 
absorption. Only four profile shape parameters were 
refined and these were the absorption coefficient //, the 
receiving slit width w,, the angular aperture of the inci- 
dent beam Soller slits A^ and the apparent crystallite 
size r^pp. All the other instrument parameters were fixed 
at their actual values. Some parameters, such as the 
equatorial divergence a, the width of the x-ray source 
and the sample width L^, have very little effect on the 
profile shapes and a wide range of values can be 
accommodated for each of these parameters without 
affecting the quality of the fit. Although the receiving 
slit length L, and x-ray source length L^ can have a sig- 
nificant effect on the profile shape, they were not 
refined because of their strong correlation with the 
Soller slit aperture A^ and the fact that the improvement 
in fit with their inclusion was minimal. Equally good 
fits are obtained by fixing L^ and A^ and allowing L, to 
refine. The equatorial divergence angle a and other 
axial divergence parameters are normally included in a 
refinement when low angle diffraction lines, preferably 
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Table 3. Comparison of actual diffractometer and sample parame- 
ters with values obtained for MgO data by FPPF refinement. In this 
refinement the radius of the diffractometer was fixed at 200 mm and 
the CuKa spectrum adopted for fitting was the one given earlier in 
Table 2. Note: NR denotes not refined 



Physical variable 


Actual value 


Refined value 


X-ray source width w^ (mm) 


0.04 


0.04 (NR) 


X-ray source length L^ (mm) 


12 


12 (NR) 


Incident SoUer shts 4 (°) 


2 


1.85 


Equatorial divergence a (°) 


0.9 


0.9 (NR) 


Axial sample length Lg (mm) 


20 


20 (NR) 


Absorption jU (cm~^) 


74 


68 


Width of receiving slit Wj. (mm) 


0.13 


0.095 


Length of receiving sht L^ (mm) 


10 


10 (NR) 




Overall Fit R^ 


3.5 % 




Expected R^ 


2.4 % 



below 20 = 25°, exist in the data set as this region of the 
diffraction pattern that is sensitive to equatorial and 
axial divergence. The inclusion of crystallite size 
broadening in most refinements is very important. In 
the present MgO refinement its removal from the 
refinement resulted in the other refined parameters tak- 
ing on physically unrealistic values and a rise in the R^ 
value from 3.5 % to 6.2 %. It was also clearly evident 
that the calculated profiles no longer fitted the tails of 
the observed profiles. 



7. Concluding Remarks 

Fundamental parameters profile fitting offers a num- 
ber of benefits as a method of profile analysis. It is 
based on a physical model of the diffractometer and its 
refined parameters should be self consistent with phys- 
ical dimensions of the diffractometer and the physical 
properties of the sample. On this basis it can therefore 
identify whether or not a diffractometer is operating at 
optimum resolution for the conditions used and provide 
a means for assessing the performance of a diffrac- 
tometer in a particular application. For example, a 
knowledge of the performance of a diffractometer oper- 
ating under asymmetric conditions can be important in 
choosing the best 20 range and slits that will maintain 
sufficient resolution to monitor line positions in stress 
analysis. As the profile shape is known, the FPPF tech- 
nique also provides greater certainty in the identifica- 
tion of weak peaks or impurity lines embedded in the 
tails of stronger line. In Rietveld analysis or quantita- 
tive analysis, FPPF allows the profile shapes across the 
whole 20 range to be fitted without any instrument 
based parameters in the refinement. The focus of the 



refinement is therefore on the diffraction effects of the 
specimen and not the instrument. Moreover, line broad- 
ening analysis is an integral feature of FPPF for both 
individual lines or as part of a Rietveld analysis and, 
correction for instrument broadening is an intrinsic part 
of the analysis making reference specimens unneces- 
sary in many circumstances. FPPF has its weaknesses. 
Up until now it has only had limited success as a 
method of accurate lattice parameter determination and 
in some diffractometer set-ups the FPPF instrument 
parameters can differ significantly from the actual val- 
ues. Although it corrects the 20 positions of lines for 
instrument line shift, including zero shift and specimen 
displacement, it has not been possible to obtain a set of 
lattice parameters for a particular specimen which are 
the same, within an uncertainty < ±0.0005 A, for every 
hkl line in the pattern [43]. In addition, the FPPF 
approach does not have a physically based model for 
incorporating mirrors and monochromators into wave- 
length distribution and the axial divergence profile 
although recent work by Masson et al. [44] suggested 
an analytical approach in which devices within the opti- 
cal path of a diffractometer can be readily incorporated 
into the theoretical formalism for describing the profile 
shape. The current physical modelling of laboratory 
based diffractometers and some of the diffraction 
processes within FPPF is adequate for many applica- 
tions, but some degree of refinement of the theoretical 
procedures underpinning the technique is still neces- 
sary before the technique can claim to accurately 
describe both the shape and position of powder diffrac- 
tion lines. 
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